Machine learning intelligent based hydromagnetic thermal transport under Soret and Dufour effects in convergent/divergent channels: a hybrid evolutionary numerical algorithm

In this research, we analyze the complex dynamics of hydro-magnetic flow and heat transport under Sorent and Dofour effects within wedge-shaped converging and diverging channels emphasizing its critical role in conventional system design, high-performance thermal equipment. We utilized artificial neural networks (ANNs) to investigation the dynamics of the problem. Our study centers on unraveling the intricacies of energy transport and entropy production arising from the pressure-driven flow of a non-Newtonian fluid within both convergent and divergent channel. The weights of ANN based fitness function ranging from − 10 to 10. To optimize the weights and biases of artificial neural networks (ANNs), employ a hybridization of advanced evolutionary optimization algorithms, specifically the artificial bee colony (ABC) optimization integrated with neural network algorithms (NNA). This approach allows us to identify and fine-tune the optimal weights within the neural network, enabling accurate prediction. We compare our results against the established different analytical and numerical methods to assess the effectiveness of our approach. The methodology undergoes a rigorous evaluation, encompassing multiple independent runs to ensure the robustness and reliability of our findings. Additionally, we conduct a comprehensive analysis that includes metrics such as mean squared error, minimum values, maximum values, average values, and standard deviation over these multiple independent runs. The minimum fitness function value is 1.32 × 10−8 computed across these multiple runs. The absolute error, between the HAM and machine learning approach addressed ranging from 3.55 × 10−7 to 1.90 × 10−8. This multifaceted evaluation ensures a thorough understanding of the performance and variability of our proposed approach, ultimately contributing to our understanding of entropy management in non-uniform channel flows, with valuable implications for diverse engineering applications.

www.nature.com/scientificreports/estimates a major research focus.Bejan 1 investigated entropy optimization to gauge uncertainty (irreversibility) in advanced engineering systems, striving to enhance functionality by minimizing entropy production through critical parameter reduction.
Various extrusion techniques and engineered systems have evolved based on thermodynamic principles.The first law of thermodynamics, conserving energy within a system without losses, guides this development.However, it does not consider energy creation.Research by Mehryan et al. 2 studied entropy behavior in magnetic third-degree fluid flow over a corrugated plate, revealing that the mean Brinkman number fosters total entropy generation.Complementing the first law, the second law of thermodynamics offers insights into entropy generation, vital for resistance management.This has led to extensive studies on entropy generation phenomena.Seyyedi et al. 3 explored entropy formation in an L-shaped enclosure with nanoparticle flow, while Riaz et al. 4 examined optimized viscoelastic nanoparticle flow in annuli with flexible walls.Turkyilmazoglu 5 studied velocity slip effects in metallic channels for optimal flow.Khan et al. 6 investigated entropy production patterns in Casson nanofluid flow caused by stretched disks, and Hayat et al. 7 analyzed numerical entropy generation assessments for the Crosser model.
Furthermore, the Jaffrey-Hamel flow in diverging/converging channels has significance in various fields.Barzegar Gerdroodbary et al. 8 conducted into the impact of thermal radiation on these channels, aiming to understand its effects on thermal profiles and configurations under different flow conditions.The foundation for studying such flows was laid over a century ago by Jeffery and Hamel, simplifying the Navier-Stokes equations and exploring thermal performance in a Newtonian fluid with nonparallel walls, known today as Jeffery-Hamel flow 9 .Yarmand et al. 10 created a hybrid nanofluid of activated carbon and graphene-EG, achieving a 6.47% increase in thermal conductivity at 40 °C with a volume fraction of 0.06%.Makinde investigated fluid flow irreversibility in a channel with variable viscosity and nonuniform temperatures 11 .Beǵ and Makinde 12 explored inherent irreversibility in nonuniform channels.Further investigations on entropy optimization in different nanofluids are found in references [13][14][15][16][17] .
Rehman et al. 18 computed entropy generation in non-Newtonian and Eyring-Powell nanofluid flows from a stretching surface.Bejan 19 derived a mathematical expression for minimizing the entropy generated in engineering systems.Due to a number of important emerging technologies, e.g., electroconductive materials processing, magnetic nozzle design, bio-inspired propulsion etc., boundary layer studies with entropy generation have stimulated extensive interest in recent years.The investigation of solutal and thermal movement in a porous medium has garnered significant attention in both theoretical and practical research, leading to notable applications in various fields, including energy storage units, geothermal systems, nuclear waste repositories, heat insulation, catalytic reactors, and drying technologies 9 .Kumar et al. 20 further extended the understanding of thermal diffusion and radiation effects on unsteady magnetohydrodynamics (MHD) flow through a porous medium.The researchers considered variable temperature and mass diffusion, and also accounted for the presence of a heat source or sink, making their investigation more comprehensive.Magnetohydrodynamics (MHD) is a specialized field that explores the interaction between magnetic fields and moving, conducting hybrid nanofluids 21,22 .
When an electrically conducting fluid is subjected to a magnetic field, it experiences a force known as the Lorentz force.This force is proportional to the fluid velocity and always opposes the flow, acting as a damping effect.Interestingly, a less widely known approach to generate a force within a flowing fluid is through the application of both an externally applied magnetic field and an externally applied electric field.This combination results in the generation of the Lorentz force, which can be achieved by arranging flush mounted electrodes and permanent magnets with alternating polarity and magnetization.The Lorentz force, when acting parallel to a flat plate, can either assist or oppose the flow.The concept of using the Lorentz force to stabilize a boundary layer flow over a flat plate can be attributed to the work of Henoch et al. 23 .Their contributions have significant implications for flow control and stability enhancement in various engineering applications.The significant role of mass and energy fluxes, arising from temperature and concentration gradients, respectively, lies in various applications such as chemical processing equipment design, crop frost damage, and fog formulation.The phenomena of Dufuor effect (diffusion-thermo) due to concentration gradients and Soret effect due to temperature gradients are studied.Analytical and numerical techniques, including the Homotopy analysis method (HAM) 24,25 , Adomian decomposition method (ADM) 26,27 and numerically 28 are utilized to solve the nonlinear governing equations.Jeffery-Hamel and spinning Disk addressed considering thermal convection considering nonlinear radiative heat transfer [29][30][31] .
In a recent year, researchers have increasingly focused their efforts on seeking solutions to complex problems in various engineering and sciences domains through the application of machine learning techniques.Traditional analytical and numerical methods have struggled to solve non-linear system of differential equations, but machine learning techniques have proven capable of providing solutions to such challenging problems with more accuracy.The bio-heat equation is an essential tool for studying heat transfer and thermoregulation in living tissues.Traditional approaches have employed deterministic solvers to evaluate the dynamics of the bio-heat equation.However, recent advancements in stochastic optimization, particularly in the context of artificial intelligence, have opened up new avenues for solving complex differential equations, including those with fractional derivatives 32,33 .Stochastic optimization techniques have found applications in various scientific domains, such as astrophysics 34 , plasma physics 35 , cell-growth modeling 36 , fluid dynamics [37][38][39][40] , and many others.Bio-inspired artificial intelligence, particularly genetic algorithms (GAs), has emerged as a powerful approach for stochastic optimization.GA is global optimization methods and have been widely utilized to solution a diverse of nonlinear problems in physiological sciences [41][42][43][44][45] .These methods offer numerous advantages, including ease of implementation, broad applicability, stability, avoidance of divergence, and high reliability.
Machine learning is an exceptionally potent tool especially artificial neural network (ANNs) renowned for its precision in solving intricate non-linear problems.The hydromagnetic thermal transport under Soret and Dufour effects in convergent/divergent channels is computed using traditional analytical and numerical approaches 46 .The research gap is thermal transport under Soren and Dofour effect in both convergent/divergent channels www.nature.com/scientificreports/through machine learning techniques is unutilized to address this phenomena.In this research we used artificial neural networks (ANNs) integrate with nature inspired evolutionary optimization algorithm artificial bee colony (ABC) hybridization with neural network algorithm (NNA) as ANN-ABC-NNA to tackle this problem and compare our finding with others traditional approaches.The statistical comparison between the ANN-ABC-NNA approach and traditional methods are computed to effectiveness of proposed techniques.

Problem formulation
The two-dimensional Carreau liquid flow between two crossing plates that are indefinitely long in the z-direction is shown in Fig. 1.The source of the flow is situated at the intake where two plates converge.Assume that there is a 2 α angle between the walls.We assume that the direction of the flow drift is radial.Moreover, the fluid veloc- ity in both the divergent ( β > 0) and convergent ( β < 0) channels is greatly impacted by the lubrication of the channel wall.In the mathematical formulation of the problem, we consider the flow to be laminar, stable, and incompressible.Assume that the flow field is subjected to an external magnetic field B 0 that is applied perpen- dicularly and has a significant impact on fluid movement.Let us consider a thermal flow under the magnetic effects in a wedge-shaped convergent/divergent channel shown in Fig. 1.The mathematical model for the pure radial physical problem ( U θ = 0 ) can be expressed as the set of following partial differential equations 46 .
Mass balance equation
Symmetries in the middle line θ is as follow boundaries: where v f , µ 0 , ρ f , k f , σ , c p , C s , σ * , k * , D B , K T , D r , designate the kinematic viscosity, dynamic viscosity, density, the thermal conductivity, electrical conductivity, heat capacitance, concentration susceptibility, Stefan-Boltzmann constant, mean absorpation coefficient, Brownian diffusion, thermal diffusion and thermophoresis diffusion.γ represents fluid-wall friction, ranging from smooth (γ = 0) to rough (γ → ∞), and δ denotes temperature slip factor in channel walls and Q the flow rate of channel in the integral form as where Q greater than zero for a divergent channel and Q less than zero for a convergent channel.Entropy (E) generation equation [1][2][3][4][5][6][7]46 : To dimensionally standardize the subsequent boundary value problem, the following transformation is used: ( where dimensionless parameters are given as here We, Re, M, Pr, R, Ec, Df, Sc, Sr, , A and m represented Weissenberg number, Reynolds number, magnetic number, Prandtle number, radiation parameter, Eckert number, Dufour Parameter, Schmidt number, Soret number, Skin friction coefficient, diffusion parameter, temperature slip and the friction wall coefficient respectively.Br (PrEc) is the Brinkman number, is the, R d is the molar gas constant, is the, and is.Nusselt number Nu, and Sherwood number Sh quantify nanofluid's flow, heat, and mass transfer rates at a boundary, representing generalized forms of physical quantities.

Solution of the problem
We using machine learning technique integrated with meta-heuristic algorithm for the solution of system of nonlinear differential equations for hydro-magnetic thermal transport under Soret and Dufour effects in convergent and divergent channels.The solution methodology is following steps.
• Neuro-computing based mathematical formulation of system of non-linear differential equation

Neuro-computing based model
In feed-forward artificial neural networks (ANNs), a prevalent strategy for approximating the solution f , θ, � and their nth order derivatives denoted as 19), ( 20) and ( 21) revolves around the incorporation of continuous mappings.These continuous mappings form the core architectural elements of the neural network, facilitating its ability to effectively replicate the desired approximation solution of f , and .The following equations have been converted to neuro-computing based model.
The ANNs based temperature and concentrations are following equations: (18) Vol.:(0123456789) are set of corresponding ANNs weights.When using these equa- tions in feed-forward artificial neural networks, the expressions of fitness function for velocity, temperature and concentration equations and associated boundary conditions are transformed.( 21)

Artificial bee colony
The artificial bee colony (ABC) is meta-heuristic optimization algorithm to find the optimum value of the function.ABC works based on comprises employed, onlooker, and scout bees.Employed bees exploit and share food source info with onlookers, while scouts seek new sources.Dances by employed bees communicate food quality; onlookers choose sources based on dance probabilities, when a source is depleted, employed bees become scouts, reflecting exploration-exploitation dynamics.In the ABC algorithm, sources represent solutions, nectar indicates fitness, and onlookers choose based on probabilities.The ABC process involves four basic phases: (a) initialization, (b) employed bee, (c) onlooker bee, (d) scout bee.
Initialization of the population.To begin with ABC initiates by generating a population solutions distributed uniformly.Every solution, denoted as x i comprises a D-dimensional vector, with D representing the number of weights associated with the Neuro-computing based model.Each x i corresponds to an individual food source within the population.The creation of each food source adheres to the following pattern: where x j min and x j max are bound of x i in j th direction.
Employed bees phase.In this phase, worker bees modify the current solution by incorporating their personal experiences and assessing the Neuro-computing based fitness function of the potential new solution.If the newly discovered food source demonstrates a higher fitness value compared to the existing one, the bee relocates to the new position, abandoning the previous location.The formula governing the update of the position for the i th candidate in the jth dimension during this phase is articulated as follows: here the term ∅ ij x ij − x kj representation to as the step size.Here,k ∈ 1, 2, . . .., Pop , and j ∈ {1, 2, . . ...D} represent two randomly selected indices and k must differ from i to ensure that the step size has a substantial impact, and ∅ ij belong to [− 1, 1].
(24) www.nature.com/scientificreports/Onlooker bees phase.In this phase, all the employed bees relay crucial fitness data regarding their improved solutions and disclose their precise positions to the onlooker bees residing within the hive.The onlooker bees, upon receiving this information, engage in a thorough analysis and make their selection of a solution based on a probability referred to as P i .This probability P i is directly correlated with the fitness of the solutions.
here fit i denoted the fitness value of the ith solution.Like a bee, the onlooker bee updates its memory stored position and assesses the fitness of a candidate, adopting the new position if it proves superior to the previous one while discarding the old position otherwise.
Scout bees phase.When a food source remains stationary for a set number of cycles, it's deemed abandoned, triggering the scout bees phase in which the abandoned source ( x i ) is replaced by a randomly chosen one from the search space in the ABC algorithm, with this cycle count known as the "limit for abandonment." where x j min and x j max are bound of x i in jth direction.

Neural networks algorithm
Neural network algorithm (NNA) 66 this innovative meta-heuristic approach blends concepts from both artificial neural networks (ANNs) and biological nervous systems.While artificial neural networks are primarily designed for predictive tasks, NNA cleverly integrates neural network principles with randomness to address complex optimization various scientific problems.Utilizing the inherent structure of neural networks, NNA demonstrates robust global optimum capabilities.Remarkably, NNA distinguishes itself from traditional meta-heuristic methods by relying exclusively on population size and stopping criteria, eliminating the need for additional parameters 66 .NNA algorithm comprises these four crucial core components: Update population.By the scenarios of NNA the population Y t = y t 1 , y t 2 , y t 3 , . . ., y t M undergoes updates via the weight matrix W t = w t 1 , w t 2 , w t 3 , . . ., w t M where w i t = w t i,1 , w t i,2 , w t i,3 , . . ., w t i,M represents the weight vector of the ith individual, and y i t = y t i,1 , y t i,2 , y t i,3 , . . ., y t i,E signifies the position of the jth individual.Notably, E denotes the count of variables.Furthermore, the generation of a new population can be mathematically articulated as follows: here M represents the population size while t corresponds to the present iteration count.The solution for the jth individual at time t is denoted as y t j , and y t new,j signifies the solution for the jth individual at the same time point, calculated with appropriate weights.Furthermore, the weight vector w t j is subject to the following formulation: Update weight matrix.The weight matrix W t is a pivotal role within NNA process of generating a novel popu- lation.The dynamics of the weight matrix W t can be refined through: where 2 represents a random value belong to [0, 1] uniform distribution and w t obj is the objective weight vector.The important point is both w t obj and the target solution x t obj share corresponding indices.To elaborate further, if x t obj matches x t v , ( v ∈ [1, M] ) at time t , then w t obj is equivalently aligned with w t v .
Bias operator.The role of the bias operator within NNA is to bolster its capacity for global exploration.A modification factor, denoted as β 1 , assumes significance in gauging the degree of bias introduced.This factor is subject to updates via: The bias operator encompasses both a bias population and a bias weight matrix, each characterized as follows: Within the bias population operator, two variables come into play-a randomly generated number M p , (29) y t j = y t j + y t new,j , j = 1, 2, 3, . . ., M. ( along with a set denoted as P. Let l = (l 1 , l 2 , l 3 , . . ., l D ) and u = (u 1 , u 2 , u 3 , . . ., u D ) represent the lower and upper limits of the variables, respectively.M p is determined as β t 1 × E signifying the ceiling value of the product of β t 1 and E. The set P consists of M p randomly selected integers from the range between [0, E] .Consequently, the bias population can be precisely defined as: here 3 represents a random value distributed uniformly within the range of [0, 1].The bias matrix also involves two variables, namely, a randomly generated number M w and a set denoted as R .The value of M w is calculated as the ceiling of β t 1 × M .In parallel, the set R comprises M w integers randomly selected from the interval [0, M].Consequently, the bias weight matrix can be accurately delineated as: where 4 , is a random number between [0, 1] subject to uniform distribution.
Transfer operator.Transfer operator is to generate a best solution toward the current optimal solution, which focuses on the local search ability of NNA.This is represented as following equation where 5 is a random number from [0, 1] uniform distribution, like the other meta-heuristic optimization algo- rithms NNA is initialized by where 6 is a random value between [0, 1] .The flow chart of the whole study is structured in Fig. 2

Results and discussion
There are difficulties in finding an exact solution when there exists complexity in the system.To meet this challenge, the numerical solutions are obtained using artificial neural networks using the hybridization of two algorithms artificial bee colony (ABC) optimization and neural network algorithm (NNA) algorithm.).As, the system solved different analytically and numerical for comparing with those obtained by proposed method.The neuro-computing based fitness function modeled for this problem is presented in Eq. ( 26) and optimized using ACB-NNA over a range of best ANNs weights.The solutions are found using inputs ranging from [0, 1] with a step size of h = 0.1 and n = 11 .The optimized weights obtained by ABC-NNA are presented in Table 1 which is used to compute numerical solutions f .Figure 3a illustrates the comparison between traditional 67 and the ANN-ABC-NNA approaches for velocity shows higher accuracy.The numerical results for convergent/divergent channel problem obtained by proposed ANN-ABC-NNA, along with their absolute errors (AE) are presented in Table 2 and Fig. 3b.The results obtained through ANN-ABC-NNA demonstrate close proximity to the traditional approach as in Table 2 and also shown in Fig. 3a.In order to thoroughly evaluate the performance of the algorithm, a comprehensive analysis of the results was performed over hundred (150) independent runs.These runs allowed an exhaustive study of the behavior of the algorithm over multiple iterations.The fitness graph and mean squared error (MSE) show in Fig. 3c,d show the evolution of fitness scores over the course of these independent runs.
The statistical analysis of the results is based on various parameters as shown in Table 3.It provides insightful metrics such as mean value average , minimum value (min) , maximum value (max) and standard deviation (S.D.) value for hundred (150) independence runs and provides a quantitative overview of the algorithm's performance over multiple iterations.This section is dedicated to the effect of various parameters, including Reynold number (Re), magnetic field parameter (M), frictional wall parameter (m), Prandtl number (Pr), thermophoresis parameter (Nt), Dufour number (Df), Brownian diffusion parameter (Nb), radiation parameter (R), Schmidt number (Sc), the Soret number (Sr), diffusion number (Δ) and the Brinkman number (Br).The subsequent discussion will delve into the physical implications of these results and will be illustrated using Figs.4, 5, 6 and 7.
Figure 4a,b show that an improved velocity profile f (ξ ) is produced by higher Reynolds numbers (Re) in both divergent and convergent channels.This occurrence emphasizes how the fluid dynamics in these channels are influenced by inertial forces.The flow pattern is significantly impacted by an increase in Reynolds number because of a rise in inertia-driven flow and channel narrowing.Pressure is raised and flow is accelerated due to the channel's narrowing and increased inertial forces.In essence, a greater Reynolds number indicates that inertial forces have a more significant effect and cause variations in the velocity distribution.The flow properties in divergent and convergent channels are shown to be significantly altered by changes in Reynolds number, while the flow rate remains constant due to the wall friction boundary conditions.This is visually represented in Fig. 4a,b.Furthermore, the connection between pressure and viscosity is made clear.A higher viscosity requires a higher pressure gradient in order to maintain a steady flow rate.As seen in Fig. 4a,b, this requirement results in higher outflow and center inflow velocities.Understanding the behavior of fluid flow in divergent and convergent channels requires an understanding of the interaction between viscosity, pressure, and flow dynamics.Let's go on to Fig. 4c,d, which illustrate yet another intriguing facet of the conversation.For both divergent and convergent channels, there is a proportional rise in velocity with higher values of (n).This suggests that the (36) y t j,P(S) = l t P(S) + u P(S) − l P(S) × 3 , S = 1, 2, 3, . . ., M P (37) velocity distribution is significantly influenced by the channel shape, which is represented by the exponent n.These graphic depictions make the relationship between channel shape and velocity clear and offer important insights into how these parameters are interdependent.The velocity increase with increasing angle in the diverging channel shown in Fig. 5a suggests that channel geometry has a major impact on fluid dynamics.Higher velocities are encouraged by the divergence, maybe as a result of the flow region expanding and enabling the fluid to accelerate.On the other hand, when the angle is decreased in the convergent channel, Fig. 5b shows a drop in velocity.This observation implies that the fluid flow is constrained by the convergent geometry, leading to reduced velocities.The fluid is squeezed in the channel created by the convergent design, which lowers the flow rate.The influence of the magnetic number (M) on the velocity profiles in divergent and convergent channels is illustrated graphically in Fig. 5c,d.The fluid close to  www.nature.com/scientificreports/ the wall accelerates in convergent channels to speeds higher than the centerline speed.Magnetic forces, which interact with the fluid and change its flow behavior, may be to blame for this acceleration.The pictures clearly illustrate how magnetic fields can obstruct flow separation in divergent channels.A smooth and stable flow is maintained by preventing flow separation, which is an important factor that can be useful in a variety of applications.Moreover, the discussion highlights the fact that in both divergent and convergent channels, higher values of the magnetic parameter (M) cause streams to concentrate near the channel center.This stream concentration is a remarkable phenomenon that is correlated with a decline in wedge flow because of increased Lorentz pressures.The resistance of nanoparticles is mostly determined by the increased Lorentz force.In both kinds of channels, the increased resistance to nanoparticles has an impact on the Carreau liquid velocity.The Prandtl number (Pr) and its impact on temperature decrease are highlighted in Fig. 6a.One important element is the reversible effect of the Prandtl number on thermal conductivity.Temperature decreases with increasing Prandtl numbers.Higher Prandtl numbers are thought to be responsible for this occurrence since they can both thin the thermal layer in divergent channels and improve heat transport.The dynamic character of heat transmission under various conditions can be seen from the link between temperature and Prandtl number.
The analysis is expanded to include Dufour numbers (Df) and how they affect temperature field solutions in Fig. 6b.The findings suggest that greater fluid temperatures in divergent channels are correlated with a rise in the Dufour number.The Dufour coefficient clarifies the intricate relationship between heat and mass transfer www.nature.com/scientificreports/by providing a concise description of the diffusion-thermo effects.It is necessary to comprehend how changes in the Dufour number impact temperature fields in order to forecast and enhance the heat transfer processes in divergent channels.The Prandtl and Dufour numbers, among other physical parameters, are important in determining temperature profiles.When working on systems involving heat transport via divergent channels, engineers and researchers can get significant insights from the reversible effects of Prandtl numbers and the subtle impact of Dufour numbers.Improved control and optimization of thermal processes in various applications are made possible by this thorough understanding of the physical characteristics.The concentration profile is shown in Fig. 7a,b, which highlight the effects of two important parameters: the Prandtl number (Pr) and the Dufour number (Df).These variables affect the system's overall behavior by playing important roles in the processes of heat and mass transport.A dimensionless quantity known as the Prandtl number (Pr) expresses the proportion of momentum diffusivity to thermal diffusivity.Momentum diffusivity predominates over thermal diffusivity when the number is larger.Figure 7a highlights a significant tendency in the context of the concentration profile: an increase in the Prandtl number (Pr) is correlated with an elevation in the concentration profile.This implies that the impact of the momentum diffusivity is stronger, resulting in a more intense distribution of concentration inside the system.Thus, knowing and adjusting the Prandtl number (Pr) might be a useful tactical tool to manage and enhance concentration profiles in the scenario under study.On the other hand, another dimensionless parameter that is essential to understanding the behavior of the system is the Dufour number (Df).It shows the proportion of mass diffusion to heat diffusion.A drop in the Dufour number (Df) is correlated with an increase in the concentration profile, as Fig. 7b illustrates.This suggests that the concentration distribution decreases when heat diffusion becomes more important than mass diffusion.Because of this, the Dufour number (Df), whose lower values encourage a more concentrated distribution, becomes an important consideration for customizing concentration profiles to satisfy particular needs.
The effect of the Soret number (Sr) on the temperature field reduction is distinctly shown in Fig. 8a.Understanding the impact of the Soret number (Sr), a dimensionless parameter that describes the thermal diffusion in a mixture, on the temperature field is essential to comprehending phenomena related to mass and heat transfer.The temperature field noticeably decreases as the Soret number (Sr) rises, as seen in the image.This finding Table 2. Comparative analysis for velocity between HAM 66 , Moradi et al. 67 , Shukla et al. 14 , Keller box 46 and ANN-ABC-NNA for parameters Re = 110, β = 3°, We = 0, n = 1, M = 0. suggests that a larger Soret number (Sr) causes thermal diffusion effects to be more prominent, which in turn causes a greater systemic temperature drop.Examining Fig. 8b, we can see some interesting findings from the assessment of particle concentration under a diverging channel.A heat gradient causes an increase in the Soret number (Sr), which is responsible for the heightened species concentration seen in the picture.A key factor in changing fluid concentration is the Soret effect, which results from the combination of concentration and temperature gradients.The concentration spike shown in Fig. 8b is explained by the increased molar mass diffusivity linked to the better Soret impact.Because of this increased diffusivity, species in the fluid are likely to be more mobile, which intensifies the concentration effect.The concentration spike is an obvious indication of how the Soret number (Sr) directly affects the system's ability to enhance fluid concentration.These results highlight how important the Soret number (Sr) is in determining temperature and species concentration fields.The correlation that exists between elevated species concentration and an increased Soret number (Sr) highlights the significance of taking Soret effects into account while comprehending and forecasting mass transport and heat behaviors in systems that have thermal gradients.This improved knowledge may have ramifications for a range of applications, including environmental research and industrial activities where exact control over temperature and species concentration is essential.For entropy generation (S G ) and Bejan number (B e ), Fig. 9a-d has been included, respectively.The fluctuation in entropy generation with respect to the channel angle factor (β) and the magnetic field (M) is depicted in Fig. 9a,b.One of the most important factors affecting the system's entropy creation is the magnetic field ( M ).An rise in the system's entropy is found in direct proportion to the strength of the magnetic field.This is explained  by the complex interaction between the thermodynamic parameters of the system and the magnetic field.It is possible that the magnetic field adds more complexity to the system, which increases the creation of entropy.In a similar vein, the entropy production is also greatly aided by the channel angle factor ( β ).An increase in system entropy is correlated with a rise in the channel angle factor.This pattern could be the result of changed flow patterns or higher fluidic resistance brought on by changes in the channel angle.Controlling the channel angle factor becomes essential to entropy management in the system.Moving on to Fig. 9c,d, which shows the Bejan number profile for the channel angle factor ( β ) and the magnetic field (M), it is important to observe that the two variables have an inverse connection with the Bejan number.The system's energy transfer efficiency and the irreversibility of its heat transfer processes are both shown by the Bejan number.As the magnetic field intensity grows, the efficiency of energy transfer diminishes, leading to a higher irreversibility of heat transfer processes.This is implied by the inverse proportionality between the Bejan number and the magnetic field ( M ) in this con- text.This highlights the significance of meticulously regulating the magnetic field to maximize the effectiveness of energy transfer.Likewise, the inverse relationship found between the Bejan number and the channel angle factor ( β ) implies that changes to the channel geometry affect the energy transfer efficiency.A higher degree of irreversibility in the heat transfer processes is indicated by a decreasing Bejan number as the channel angle factor rises.This emphasizes how crucial it is to have a thorough grasp of channel geometry impacts in order to maximize the efficiency of energy transfer.

Conclusion
In conclusion, this research delved into the intricate dynamics of incompressible viscous fluid of hydro-magnetic flow and heat transport within convergent and divergent channels, shedding light on its critical significance in conventional system design, high-performance thermal equipment, and geothermal energy applications.Leveraging the power of machine learning, specifically artificial neural networks (ANNs), we undertook a comprehensive computational investigation.Our study focus was on unraveling the complexities of energy transport and entropy production resulting from the pressure-driven flow of a non-Newtonian fluid in these channel geometries.To enhance the performance of neuro-computing based fitness function, we adopted a hybridization approach of advanced evolutionary optimization algorithms, namely artificial bee colony (ABC) optimization and neural network algorithms (NNA).This allowed us to fine-tune the neural network's optimum weights and biases, ultimately leading to accurate predictions of dynamics.
• Comparing ANN-ABC-NNA results with established analytical and numerical methods underscored the efficacy of our approach in addressing this challenging problem.• The absolute error between the HAM and ANN approach ranging from 1.90 × 10 −8 to 3.55 × 10 −7 with optimized weights and biases through ABC-NNA and also improved results as compare to Keller box method.• Our rigorous methodology evaluation multiple (150) independent runs, provided robust and reliable find- ings.The statistical analysis of ANN-ABC-NNA algorithm we employed various metrics, including mean squared error, minimum and maximum values, average values, and standard deviation, to comprehensively assess our approach's performance and variability over 150 independence runs.• The minimum, maximum, average and standard deviation values over these independence runs are from 2.05 × 10 −9 to 4.39 × 10 −4 .• This research advances our understanding of entropy management in nonuniform channel flows, particularly when dealing with nano-materials an understanding that holds significant implications for a wide array of engineering applications.• The synergy of machine learning techniques and advanced optimization algorithms offers a promising avenue for tackling complex fluid dynamics problems, opening doors to more accurate and efficient solutions in engineering and related fields.• The results indicate that higher fluid temperatures in divergent channels are correlated with an increase in the Dufour number.• The channel narrows and inertial forces increase pressure and accelerate flow.
• The velocity of both convergent and divergent channels decreases with increasing magnetic field strength.
• There is a correlation between an increase in the concentration profile and a decrease in the Dufour number (Df).• As the Soret number (Sr) increases, the temperature field clearly drops.

Figure 2 .
Figure 2. Flowchart of hybrid artificial bee colony and neural network algorithm.
• ANN based fitness function • Hybrid meta-heuristic techniques ABC-NNA to optimize the fitness function for the best weights and biases of ANN within [− 10, 10].